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A toroidal trap combined with external time-dependent electric field can be used for implementing 
different dynamical regimes of matter waves, fn particular, we show that dynamical and stochastic 
acceleration, localization and implementation of the Kapitza pendulum can be originated by means 
of proper choice of the external force. 



PACS numbers: 03.75.Lm, 03.75.Kk, 03.75.Ss 



I. INTRODUCTION 



Exploring different geometries of potentials trapping 
cold condensed atoms is of both fundamental and practi- 
cal importance. Toroidal traps play a special role allow- 
ing for "infinite" atomic trajectories and for realization 
of quasi-one-dimensional (quasi-lD) regimes. These ad- 
vantages are relevant for designing highlyprecise sensors 
based on matter wave interferometry [l|, Q as well as for 
accurate study of such phenomena as superfluid currents, 
stability of sound waves, solitons and vortices in Bose- 
Einstein condensates (BEC's) [1,0] • Traps with circular 
geometry are also believed to be conceptually important 
for implementation of the main ideas of the accelerator 
physics at ultra-low temperatures [31 and, in particular, 
for acceleration of ultracold atoms In this last con- 
text existence of well localized wave packets, and thus 
attenuation of the dispersion, the latter being the intrin- 
sic property of a quantum systems, is of primary impor- 
tance. In the first experimental studies [2| it was shown 
that the dispersive spreading out [l[ can be compensated 
by using betatron resonances in a storage ring. An al- 
ternative way of contra-balancing dispersion is also well 
known - it is nonlinearity, leading in quasi- ID regime to 
existence of bright and dark matter solitons (see e.g. [1,0] 
and [1,S^ 
explored 



1C | , respectively) . This issue has already been 
11[ from the point of view of acceleration of 



matter waves in a toroidal trap with help of a modulated 
optical lattice, which is known to be an efficient tool for 
acceleration of matter waves [ijj • 

In this paper we propose two alternative ways of ac- 
celerating matter wave solitons - either by time varying 
or by stochastic external electric field. These new ways 
of soliton acceleration are especially relevant in view of 
radiative losses [l3| and distortions of solitons mov- 
ing in optical lattices (the effects acquiring significance 
for long trajectories). At the same time, it turns out 
that the toroidal geometry of a trap confining a BEC al- 
lows one to realize a number of other dynamical regimes, 
like dynamical localization of solitons and solitonic im- 
plementation of the celebrated Kapitza pendulum. The- 



oretical description of all mentioned phenomena can be 
put witching the unique framework, based on the pertur- 
bation theory for solitons, what is done in the present pa- 
per. More specifically, in Sec. [IT] we formulate the model 
and the main physical constraints determining its valid- 
ity. In sections IIIII and IIVI we describe how by applying 
external time-dependent electric field matter solitons can 
be accelerated in the usual sense and in the sense of the 
time increase of the velocity variance (the stochastic ac- 
celeration) , respectively. In Sec. [V] we describe localized 
states of the matter in circular trap subject to external 
field, and in Sec. I VII we show that a matter soliton af- 
fected by rapidly varying force represents an example of 



the Kapitza pendulum 14]. Summary and discussion of 



the results are given in the Conclusion. 



II. SCALING AND THE EVOLUTION 
EQUATION 

We assume that a BEC is loaded in a circular trap, 
which in cylindrical coordinates r = (p, ip, z) is described 
by V = V c (p) + muj 2 z 2 /2, where cu z is the frequency of 
the magnetic trap in the z-direction, V c {p) is the poten- 
tial in the radial direction, forming the trap circular in 
the (x, y)-plane, and m is the mass of an atom. We also 
suppose that the BEC is subject to external electric field 
with amplitude Eq, which produces an additional poten- 
tial V ex t — —a'Efi/A, where a' is the polarizability of the 
atoms (see e.g. |l5[). If the amplitude Eq or direction of 
the field vary along some direction, say, along the x-axis, 
smoothly on the scale of the trap radius R, the potential 
energy V ex t van be expanded in the Taylor series and, af- 
ter neglecting nonessential constant, be rewritten in the 
form V ext = —ax, where a = (a' /4)d(E ) 2 /dx\ x= o and 
we restricted the consideration only to the first term of 
the expansion. In order to realize one-dimensional geom- 
etry we require torus radius to be much larger than the 
core radius r c , what allows us to define a small param- 
eter e = r c /R < 1. In order to introduce quantitative 
characteristics, we consider the normalized ground state 
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cj) of the eigenvalue problem 

— p— <f> + V c {p)<l> = e r (l>, / (^Vp=1(1) 

2m p dp dp J Q 

and define R x = J Q °° <j> 2 p 2 dp, R 2 = (J °° (f) 2 dp/p) , 
and A = (J Q °° (ffipdp) 1 ^ 2 . In the case at hand A ~ 
y/Rr2 ~ e 1 / 2 ^ and thus A « fii ~ i? 2 ~ -R. 

In the present paper we are interested in the dynam- 
ics of matter waves which spatial extension is much less 
than the trap perimeter, allowing treat them similarly to 
the matter solitons in an infinite one-dimensional trap. 
This in particular the case where the spatial size of the 
BEC excitations along the trap are of order of A, which 
is the well defined parameter and thus convenient for for- 
mulation the constraints of the theory. Indeed, now we 
can estimate the kinetic energy of the longitudinal ex- 
citations as eii = h 2 /(2mX 2 ) and require it to be much 
less than the kinetic energy of the transverse excitations, 
e r ~ h 2 / (2mr 2 ) (for the sake of simplicity here we assume 
that the size of the trap in z-direction is of order of the 
core radius: a z — ^Jhjmus z ~ r c ). Adding the require- 
ment for the energy of the two-body interactions, which 
is estimated as \g\n (where g = 4irh 2 a s /m, a s is the scat- 
tering length, n ~ N/V is a mean density, N is the total 
number of atoms and V is the effective volume occupied 
by the atoms and estimated as V ~ irXr c a z ), to be of or- 
der of en and to be much less than e r (or more precisely, 
requiring \g\n/e r ~ e), we can neglect in the leading or- 
der the transitions between the transverse energy lev- 
els and employ the multiple scale expansion |6l.[Toj 
for description of the quasi-one-dimensional evolution of 
the BEC. We also notice that subject to the assumptions 
introduced, one has the estimate N ~ e 3 / 2 i?/(8|a s |). 

In order to get an insight on practical numbers, let us 
consider 7 Li atoms (a s = —2 nm) in a trap with R = 
100 /im, r c = 5 pm and a z = 10 pm. Then e = 0.05, the 
characteristic size of solitonic excitations is A » 22 pm 
and the number of particles is estimated as N ss 140. We 
emphasize, that these estimates indicate only an order of 
the parameters. Thus, for example, a condensate of 10 2 4- 
10 3 lithium atoms satisfy the conditions of the theory. 

We will be interested in managing soliton dynamics 
by means of weak (i.e not destroying solitons) electro- 
magnetic field varying in time. Respectively, we con- 
sider a time-dependent and characterized by the esti- 
mate a ~ h 2 / (mRiX 2 ). Then, starting with the Gross- 
Pitaevskii equation, in which the external potential in 
cylindrical coordinates has the form V ex t = —apcos(ip) « 
— aR\ cos(y), and using the multiple-scale expansion one 
ensures that the BEC macroscopic wave function in the 
leading order allows factorization 

* = 7r- 1 / 4 a ; 1 / 2 e- l ^" + ^) t / 2 e^ z2 / 2a ^(r)^(i, <p), (2) 

where cu r = 2e r /fi and ip(t,ip) solves the nonlinear 
Schrodinger equation, which we write in terms of A = 



\J \g\m I \/2/kTi 2 a z Tp 1 £ = R2ip/X, and r = Ht/mX 2 
8 A 1 r) 2 A 

i &F = ~2W~ cos(kC)/(t)a + a|A|2 ^ ' (3) 

Here a =signa s , /(r) = mRiX 2 a(t) /h 2 and K — X/R2 ~ 
y/e. We choose the scaling in such a way that all terms 
in {3| are of the unity order, and in particular A — 0(1). 
This can be done, taken into account the normalization 

/V| 2 dC = 2V2^^, (4) 

L = 2tt/k, which follows from the normalization condi- 
tion for the order parameter J \^>\ 2 d 3 r — N, and consid- 
ering TV ~ a z /\a s \, what is of order of 10 3 , in a typical 
experimental setting. 

Eq. (J3j) is subject to periodic boundary conditions 
A(C,t)=A(C + L,t). 



III. ACCELERATION OF BRIGHT MATTER 
SOLITONS BY TIME-DEPENDENT EXTERNAL 
FORCE 

First we consider the acceleration, 7, which can be 
achieved due to the potential V e properly dependent on 
time. An order of magnitude of 7 can be estimated by 
taking into account that Eq. ([3]) makes sense provided 
that all terms are of the unity order. In the physical 
units this gives 7 ~ ?i 2 /(m 2 A 3 ). Then, recalling the 
above example of the lithium condensate we estimates 
7 ~ 7mm/s 2 , what is of order of the acceleration an- 
nounced in [ll| ■ It however does not provide the best es- 
timate in our case, because it is based on the ID model, 
while lowering dimensionality imposes constrains on the 
atomic density and consequently on the amplitude of the 
applied force. 

To describe the physics of the phenomenon we consider 
a BEC with a negative scattering length (a = —1). Then 
a "bright soliton" solution of Eq. ((3|) at /(r) = (or more 
precisely a periodic solution mimicking a bright soliton 
in an infinite ID system) which moves with a constant 
velocity v n can be written down as follows 0] 

A s = e-^W+^k+^C^) dn (ry(fc)(£ - VnT ), k) . (5) 

Here 6n(x,k) is the Jacobi elliptic function [161 ], k is 
the elliptic modulus parameterizing the solution. The 
frequency and the amplitude are given by: u{k) — 
(k 2 /2 - l)if{k) and rj(k) = 2K(k)/L \K{k) is the com- 
plete elliptic integral of the first kind] . The velocity of 
the soliton is quantized v n = 2irn/ L with n being integer. 

To ensure that the solution A s satisfies the scaling 
relations imposed above, we notice that the size of 
the soliton can be estimated as ir/K(k) and its small- 
ness implies that k is close to unity. In that case we 
obtain the estimates 1 — k 2 ~ 16 exp(— 27r/ v / e) and 
dn (r)(k)(C - v n r), k) w 1/ cosh(r?(fc)(C - v n r)). 
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In the limit k — > 1 quantization of the velocity does not 
play any significant role. We verified this numerically. 
For example, for L = 10 deviation of the initial velocity 
from the quantized one produces appreciable effect on 
dynamics during intervals r < 100 only if k < 0.99. 

When external force is applied, /(t) ^ 0, the velocity is 
not preserved any more, what manifests itself in evolution 
of the momentum P = (l/2i) j£ {AqA - AAq) (here A 
stands for complex conjugation of A) according to the 
law: 



dP_ 

dT 



-fir) 



cos(«C)MV 



(6) 



The external field, however, does not affect the norm: 
N = J Q L |A| 2 dC=const. It follows from (JSJ that in the 
adiabatic approximation the solution of the perturbed 
equation §S§ can be searched in the form 



.4 : 



-iWQ+V(T)'/2)r+iV(TX v fa^_ X ( T y } Q (7) 



where V(r) = dX(T)/dT is the time-dependent velocity 
of the soliton and X(t) is the coordinate of the soliton 
center. Substituting ([7]) in ([6]) and taking into account 
the parity of the functions in the integrand as well as the 
fact that all of them are periodic with the same period 
L, we obtain the equation for the soliton coordinate 



Here 



C{k) 



d 2 X 
~d^ 



K{k) 
2nE(k) 



-KC(k)f(r) sm( K X). 



(8) 



cos(6»)dn 2 



K(k) i 



e,k)de (9) 



and it is taken into account that N — 2r]E(k), where 
E(k) is the complete elliptic integral of the second kind. 

Depending on the choice of the function /(r), Eq. (j8]) 
describes different dynamical regimes. Now we are inter- 
ested in acceleration which occurs during the rotational 
movement of the soliton in the trap (i.e. X is a growing 
function) . We illustrate this acceleration using an exam- 
ple of the simplest step-like dependence /(r). To this 
end we assume that initially the soliton is centered at 
X(0) = and require /(r) to be a constant fo for time 
intervals such that the soliton coordinates X(t) € I p and 
to be zero for X(t) £ I p where the intervals I p are given 
by I P = [ip + \)L, (p + 1)L] with p = 0, 1, .... Then, as 
it is clear, the acceleration of the soliton, which is given 
by the right hand side of JHJ, is positive for all times. 
The above requirement introduces natural splitting of 
the temporal axis in the set of intervals T/ = [77,77+1] 
(I = 0, 1, ...), with r = 0, such that /(r) = for r G T 2p 
and /(t) = f for t E T 2p+1 (here X{n) = IL/2). Thus, 
our task is to find 77 . This can be done by taking into ac- 
count that during each of the "odd" intervals T 2p +i Eq. 
([5]) describes conservative nonlinear oscillator, the solu- 
tion for which is well known. During "even" intervals T 2p 
the motion is free (with a constant velocity) what means 
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FIG. 1: (Color online) The soliton velocity vs time (panels 

a, b) for the parameters k = 0.99999, L = 10.0, fo = 0.3, 
n = 0.43 (panel a, the non-quantized velocity), n = 1.0 (panel 

b, the quantized velocity), and the forms of the soliton (panels 

c, d,e) in the instants of time indicated in the panel a. In 
the panels a,b solid and dashed lines represent the velocity 
numerically computed from Eq.Q, and Eq.©, respectively. 



that the time T 2p necessary for the soliton to cross an 
interval [pL, (p + ^)L] is 



T 2p — r 2p+ i — r 2p — L/(2v 2p ), 



(10) 



where v 2p is the velocity in the point pL. During the 
time interval T2 P +i the soliton has to cross the interval 
I p . From this condition we obtain: 



T 2 p+1 — T 2p+2 — T 2p+ l — 



V2K(^/2E /(H 2p+1 +E a 



(11) 

where H 2p+ i = u| p +i/2 + E is the energy of the soli- 
ton in the point (p + 1/2)L, Eq = C(fc)/o, and v 2p+ i = 
v 2p is the soliton velocity in the same point. At the 
end of the interval T 2p+ \ the soliton velocity is given 
by i>2( p +i) — \f^{H2 P +i + Eq). Thus one computes 
that after p rotations the soliton acquires the velocity 
v 2 (p+i), which can be obtained from the recurrent rela- 
tion: i> 2 (p+i) = Jv% p + 4C(fc)/ . 

In FigfT] a,b we compare the solution, obtained from 
the perturbation theory, Eq.([8]), with numerical simula- 
tion of Eq.© for fo = 0.3. Nevertheless during the nu- 
merical simulation we used the values for T 2p and T 2p+ i 
fEqs. pU|) - pT|l ) obtained for the case of adiabatic approx- 
imation. It follows form the results presented that the 
dashed and solid lines perfectly match until r « 50.0. At 
larger times appreciable discrepancy appears. It occurs 
due to failure of the adiabatic approximation and can be 
removed by introducing temporal corrections to T2 P and 
T 2p +\. This naturally leads to an optimization problem, 
which requires numerical approach and goes beyond the 
scope of the present work. Finally we notice, that for the 
above example of 7 Li condensate the obtained accelera- 
tion is 0.36mm/s 2 . 

Comparison of the panels a and b in FigfT] shows that 
for k « 1 quantization of the velocity is not important, 
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what is also confirmed by the evolution of the solitonic 
forms depicted in the panels Fig[T]c-e. 



IV. STOCHASTIC ACCELERATION OF 
MATTER SOLITONS 

Now we concentrate on another dynamical regime - on 
the stochastic acceleration - where increase of the veloc- 
ity of a matter soliton in a toroidal trap is achieved by 
applying a fluctuating external field. To this end hold- 
ing all conditions of the applicability of the model ([3]), 
we consider the case of a stochastic force /(r), which 
is delta-correlated Gaussian process with characteristics 
(/(t)) = and (/(r)/(r')) = D5{t - r') (the angular 
brackets stand for the stochastic averaging and D is the 
dispersion) . Now the dynamics can be described in terms 
of the distribution function 

P(V,9,t) = (*($ - *(t))*(V - V(r))), (12) 

where $(r) = kX is the angular coordinate of the soli- 
ton, <&(t) and V(r) with explicit time dependence stand 
for the soliton coordinates obtained from the dynamical 
equations while $ and V are considered as independent 
variables. The distribution function solves the Fokker- 
Planck equation, which is obtained by the standard pro- 
cedure (see e.g. |17j): 



dV 



-V 



dV 



Z?sin 2 ($) 



d 2 v 

dV 2 



(13) 



Here D = K 4 C 2 (k)D is the diffusion coefficient. Due to 
the circular geometry of the trap Eq. (fl~3|) is considered 
on the interval — it < <!> < tt with the periodic boundary 
conditions V(V, $ — it, t) = V{V : $ + tt, t) with respect 
to $ and zero boundary conditions with respect V: V — > 
as V — > ±oo. The normalization condition for the 
probability density reads: dV d^V = 1. 

Multiplying Eq. fT^)) by V and $ and integrating over 
V and <!> one readily obtains that the averaged velocity 
and angular position of the soliton are constants, which 
for the sake of simplicity will be considered zeros, i.e. 
(V) = and ($) = 0. Next, multiplying (13]) by V 2 , $ 2 
and V$ and performing the integration one obtains the 
equations of the second momenta. They are not closed 
and can be written down as follows: 



— (V 2 ) = 2L>(sin 2 $) 
|-($ 2 ) = 2<V$), 



(14) 
(15) 



d f°° 

_(y$) = -27r / P(w, V, r)V 2 dV + (V 2 ) . (16) 



Eq. (|14p means that the average square velocity is grow- 
ing with time, i.e. the soliton undergoes the stochastic 
acceleration. The law of the growth of the velocity in- 
variance deviates form the linear, as it would happen for 
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FIG. 2: (Color online) The mean square velocity (panel a) and 
the stochastic acceleration 7 (panel b) of the soliton vs time 
for different values of the dispersion, obtained by numerical 
integration of Eq. (|13[1 with parameters L = 10.0 and k — 
0.99999. In panel b dashed lines depict the approximation of 
numerical data by the law 7 oc r~ 1//2 . All axes in panels a,b 
are represented in logarithmic scale. 



the Brownian diffusion in the momentum space, what 
happens because the diffusion coefficient in the Fokker- 
Planck equation (|13[) is not a constant, but depends on 
the angular variable. However, due to the diffusion one 
can expect that the phase distribution will tend to ho- 
mogeneous, i.e. that V — > l/(27r) as r — * 00. In this 
formal limit one obtains that (V"$) — * 0, (sin 2 $) — > 1/2 
and hence (V^ 2 ) — ► Dt. In other words, the system (fT4|) - 
(Unj), describes random walk which in the limit of large 
time, approaches the Brownian diffusion in the velocity 
space. In that limit the stochastic acceleration, which 
can be defined as 7 = dyj (V 2 )/dT, would tend to zero 
according to the law 7 oc t~ x / 2 . 

In order to check the above predictions and reveal other 
features of the stochastic dynamics of a soliton in a ring 
trap we solved numerically Eq. (|13p subject to the ini- 
tial condition V(V,$,0) = (6($)5(V)}. The results is 
summarized in Fig. [5] In the panel a one observes the 
predicted monotonic growth of the mean velocity with 
time, which slightly different from the linear law. In the 
panel b one can see that the stochastic acceleration 7 
is a monotonically decreasing function, which at suffi- 
ciently large times tends to zero. In particular, at r > 15 
the decreasing of the acceleration with time can be well 
approximated by the predicted law 7 oc r^ 1 / 2 , as it is 
shown by dashed curves in the panel b of Fig. [2] (it was 
verified that in at the same times (sin 2 $) ps 1/2, what 
is in agreement with the analytical predictions). Also 
Fig. [2] b shows that the stochastic acceleration is larger 
for larger D. The physical explanation of this last fact 
is simple: the acceleration is generated by the stochastic 
force, whose intensity is determined by the dispersion D. 
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V. LOCALIZATION OF MATTER INDUCED BY 
THE EXTERNAL FIELD 

Let us now turn to localized states of a matter in a 
toroidal trap and concentrate on the states generated by 
the constant external electric field, i.e. by /(r) = /o- 
Respectively, we look for stationary solutions of Eq. (|3|) 
in the form A — e~ tUT A{C,) and obtain for A(() the equa- 
tion: 



I d 2 A 



focos(n()A + a\A\ 2 A = ujA (17) 



which is subject to periodic boundary conditions 

Several lowest branches of the numerically obtained 
solutions of Eq. (fT7|) are shown in Fig. [3] The lowest 
branch approaches zero at the frequency loq rs —0.143 
(it is interesting to mention that this frequency coincides 
with the lowest gap edge of the spectrum of the Mathicu 
equation (|17[) considered on the whole axis), where the 
amplitude of the nonlinear periodic mode is small and 
it transforms into the linear periodic Bloch mode at the 
lowest gap edge. Such a behavior of the branch is sim- 
ilar to that of the st rong ly localized modes in a BECs 
in the optical lattice [18|. The lowest mode A is local- 
ized in the vicinity of <p — (FigJJl b), i.e. around the 
minimum of the effective potential and that is why such 
a mode is stable and can exists even in the linear case, 
where the two-body interactions are negligible (here it 
is important that we are dealing with periodic boundary 
conditions). The modes of the upper branches - B and 
D (their typical examples are shown in FigJ3]b) - bifur- 
cate at lu* rs —0.345. They are localized at <p = ±7r, i.e. 
near the points where the potential has maxima. As it 
is clear, this is pure nonlinear effect and occurs due to 
delicate balance between the attractive interactions and 
repulsive forces of the external potential. Such balance 
can easily be destroyed even by an infinitesimal pertur- 
bation, what allows us to expect instability of the modes. 
The mode C represents two local maxima of the atomic 
distributions at ip = 0, it. Similar to the modes B and D 
we one can expect it to be unstable, what physically can 
be explained by existence of local atomic maxima at the 
maxima of the potential. By direct numerical solution of 
Eq. ^ (more specifically by perturbing the mode pro- 
files by the factor 1 + 0.1 cos (21 £) and computing the 
dynamics until r = 1000) we have verified that, indeed, 
only the mode A on Fig. [3] is dynamically stable, while 
the modes B, C, and D are unstable. 



VI. MATTER SOLITON AS A KAPITZA 
PENDULUM 

As the final example of nontrivial dynamics of a mat- 
ter soliton in a toroidal trap we consider dynamical lo- 
calization induced by a rapidly oscillating force f(r) = 
fa [v + cos(Or)]. In this case the solitonic motion mimics 
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FIG. 3: (Color online) The number of bosons N (for the ex- 
ample of the lithium condensate described in the text) us 
frequency u) (panel a) and shapes of the localized modes at 
iv — —0.4 (panel b) for the case where L — 10.0, /o = 0.3, 
a = -1. 



the famous Kapitza pendulum, which acquires an addi- 
tional stable point due rapid oscillation of the pivot [l4[ . 
Assuming that the physical conditions of the validity of 
the quasi- ID approximation holds and that the fre- 
quency ft is large enough, i.e. tt 2 ^> K 2 C(k)fo, one can 
perform the standard analysis (see e.g. [Hj]), i.e. look for 
a solution of in a form X(t) + £(t) where £ is small, 
|£| -C \X\, and rapidly varying, and provide averaging 
over rapid oscillations. Then one arrives at the equation 
d 2 X/dr 2 = —dU/dX with the effective potential 



U = -C(k)f 



vcos(kX) + K 2 C{k)^ cos(2kA) 



.(18) 



If the condition k 2 C (k) fo / \2£l 2 v) > 1 is met, the ef- 
fective potential U possesses two stable points: X = 
($ = 0) and X = L/2 ($ = n). So, it opens the pos- 
sibility for the new type of soliton moving around the 
new stable point. Two typical trajectories of the soli- 
ton, obtained by numerical integration of Eq. ([3]), are 
presented in Figf?] One of the trajectories displays oscil- 
lations around the new equilibrium point, while the other 
one shows the large oscillations around the equilibrium 
point $ = 2tt started with the same initial data but in the 
case where $ = tt is not an equilibrium any more. The 
amplitude of large oscillations decay with time because 
of energy losses of the soliton in the nonconservative sys- 
tem. 



VII. CONCLUSIONS 

In the present paper we have shown that dynamics of 
a matter soliton in a toroidal trap, well reproducing one- 
dimensional geometry, can be very efficiently governed by 
time varying external electric field. In particular, such 
regimes like dynamical acceleration, stochastic acceler- 
ation, localization and implementation of the Kapitza 
pendulum can be realized by proper choices of the time 
dependence of the external force. 

Experimental detection of the acceleration can be im- 
plemented either by direct imaging of the atomic cloud, 
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FIG. 4: (Color online) The angular coordinate of the soliton 
center vs time for the soliton motion affected by the rapidly 
oscillating external force, obtained numerically from Eq.([3]) 
with parameters L — 10.0, n = 0.01, a = —1, fo = 0.15, 
n = 2.0, and k = 0.99999. 



which is well localized in space and has well specified tra- 
jectory, or by measurement of the atomic distribution in 
the momentum space displaying shift of the maximum 
towards higher kinetic energies. Alternatively, one can 
study the evolution of the atomic cloud releasing from 
the trap (by switching the trap off) after some period of 
accelerating motion. The respective dynamics will be a 
spreading out cloud whose center of mass is moving with 
the acquired velocity. 

The obtained results were based on the one- 
dimensional model, although deduced using the multiple- 
scale method and thus mathematically controllable. This 



means that a number of problem, are still left open. One 
of them is the limitation on the soliton velocity, and thus 
acceleration, introduced by lowering the space dimension, 
which appears when the solitonic kinetic energy becomes 
comparable with the transverse kinetic energy. Another 
limitation on the soliton acceleration emerges from the 
fact of the velocity quantization, when the radius of the 
ring trap is not large enough. We also left for further 
studies the diversity of localized stationary atomic distri- 
butions supported by the external filed, indicating only 
the lowest modes. We thus believe that the richness of 
the phenomena which can be observed by simple com- 
bination of the trap geometry and varying external field 
will stimulate new experimental and theoretical studies. 



Acknowledgments 

The authors are indebted to H. Michinel for providing 
additional data on Ref . [ll| . YVB was supported by the 
FCT grant SFRH/PD/20292/2004. VVK was supported 
by the Secretaria de Stado de Universidades e Investi- 
gaci?n (Spain) under the grant SAB2005-0195. The work 
was supported by the FCT (Portugal) and European pro- 
gram FEDER under the grant POCI/FIS/56237/2004. 



[1] S. Gupta, K. W. Murch, K. L. Moore, T. P. Purdy, and 

D. M. Stamper-Kurn Phys. Rev. Lett. 95, 143201 (2005) 
[2] K. W. Murch, K. L. Moore, S. Gupta, and D. M. 

Stamper-Kurn, Phys. Rev. Lett. 96, 013202 (2006) 
[3] J. Javanainen, S. M. Paik, and S. Mi Yoo, Phys. Rev. A 

58, 580 (1998); L. Salasnich, A. Parola, and L. Reatto, 

Phys. Rev. A 59, 2990 (1999); A. B. Bhattacherjee, E. 

Courtade, and A. Arimondo, J. Phys. B, 37, 4397 (2004). 
[4] J. Brand and W. P. Reinhardt, J. Phys. B: At. Mol. Opt. 

Phys. 34, L113 (2001) 
[5] O. Dutta, M. Jaaskelainen, and P. Meystre, Phys. Rev. 

A 74, 023609 (2006). 
[6] V. M. Perez-Garcfa, H. Michinel and H. Herrero, Phys. 

Rev. A 57, 3837 (1998). 
[7] F. Kh. Abdullaev, et. al. Int. J. Mod. Phys. B 19, 3415 

(2005) 

[8] T. Tsuzuki, J. Low Temp. Phys. 4, 441 (1971). 

[9] Th. Busch and J. R. Anglin, Phys. Rev. Lett. 84, 2298 
(2000); A. Muryshev, G.V. Shlyapnikov, W. Ertmer, 
K. Sengstock, and M. Lewenstein, Phys. Rev. Lett. 89, 
110401 (2002) 



[10] V. A. Brazhnyi and V. V. Konotop, Phys. Rev. A 68, 
043613 (2003). 

[11] A.V. Carpentier and H. Michinel, |cond-mat/0610047| 
[12] V. A. Brazhnyi, V. V. Konotop, and V. Kuzmiak, Phys. 

Rev. A 70, 043604 (2004) 
[13] A. V. Yulin, D. V. Skryabin, and P. St. J. Russell, Phys. 

Rev. Lett. 91, 260402 (2003) 
[14] see e.g. L. D. Landau and E. M. Lifshitz, Mechanics 

(Pergamon Press, New York, 1976) 
[15] C.J. Pethick and H. Smith, Bose-Einstein condensation 

in dilute gases (Cambridge, University Press, 2001) 
[16] D. K. Lawden: Elliptic Functions and applications 

(Springer- Verlag, New York Inc., 1989) 
[17] see e.g. V. V. Konotop and L. Vazquez Nonlinear random 

waves (World Sci., Singapore, 1994) 
[18] A. Trombettoni and A. Smerzi, Phys. Rev. Lett. 86, 2353 

(2001); F. K. Abdullaev, B. B. Baizakov, S. A. Dar- 

manyan, V. V. Konotop, and M. Salerno, Phys. Rev. A 

64, 043606 (2001). 



